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SUMMARY 

A parametric analysis has been conducted to determine the interactions of 
several design variables in obtaining maximum turbine efficiency. Alkali-metal 
vapor , axial-flow turbines were investigated within a framework of several 
generalizing assumptions. Efficiency reductions due to the presence of mois- 
ture have not been included. Trends rather than magnitude were emphasized. 

Turbine tapering (i. e. , changing the ratio of last- to first-stage blade 
mean-section diameter) was an effective means for attaining maximum overall 
static efficiency with a minimum number of stages. The turbine taper ratios 
for maximum efficiency solutions decreased with increasing working fluid molec- 
ular weights. At an inlet temperature of 2500° R, these solutions were char- 
acterized by 12 stages for sodium with a taper ratio of at least 1.5, 7 stages 
for potassium with a taper ratio of roughly 1. 3, 5 stages for rubidium with a 
taper ratio of 0.985 (nearly constant mean diameter), and 3 stages for cesium 
with a taper ratio of 0.95 (nearly constant tip diameter). 

Changes in turbine taper ratio and/or changes in stator exit angles were 
shown to be a possible means for providing wide margins of selection in design 
rotational speed, while still maintaining a high level of efficiency. The ad- 
vantages of turbine tapering in reducing the number of stages required for 
nearly maximum efficiency were shown to decrease with decreasing inlet tempera- 
ture. The obtainable maximum efficiency decreased with decreasing ratios of 
outlet to inlet temperature. Hence, selection of a low ratio of outlet to 
inlet temperature may penalize turbine and system performance. 

Over a generator power output level from 1 to 10 megawatts, it was found 
that turbine diameter and rotational speed could be scaled directly within an 
error of ±2 percent, independent of any prescribed turbine taper ratio and 
number of stages. Two working fluid states were investigated. Ramifications 
of the choice of disk taper ratio limit, last-stage hub-to-tip radius ratio, 
blade-root stress correction factor, and small changes in assumed stress level 
were explored. 



IimODUCTION 


Turboelectric propulsion is currently being considered for future space 
missions. In this type of system, nuclear energy is converted to mechanical 
energy in a high-temperature Rankine cycle using an alkali metal as the working 
fluid. The turbine in this cycle drives an alternator, which produces the 
electrical energy that is used for propulsion. One requirement of such sys- 
tems is low specific powerplant weight. 

Although the turbine is comparatively light (see ref. l), its performance 
has a direct effect on the size of the radiator and, therefore, on the weight 
of the system. The gross effects of turbine inlet temperature and efficiency 
on required radiator area are shown in reference 2. For example, a 200° in- 
crease in inlet temperature reduces radiator area by about 30 percent, while a 
10-percent decrease in efficiency increases the required radiator area by about 
15 percent. This reference also contains an estimate of the number of stages 
required to achieve high efficiency with a stage speed- work parameter of 0.5 
for working fluids of rubidium, potassium, and sodium. Reference 3 presents a 
more detailed look at turbine characteristics at the 1-megawatt power level. 

The efficiency correlation of reference 4 is used to show the minimum number 
of stages required to obtain a level of at least 0.8 in aerodynamic efficiency. 
It is also shown in reference 3 that the use of increasing rather than constant 
blade mean-section diameters from the first to the last stage results in fewer 
required stages. 

The purpose of this investigation was to take a more detailed look at the 
turbine as one of the Rankine cycle system components. In general, the turbine 
designer may select among many geometric parameters to best fit the thermo- 
dynamic requirements for each particular application. On the other hand, the 
systems engineer must specify the thermodynamic constraints to best fit the 
overall powerplant requirements. Both of these problems have been investi- 
gated; emphasis was placed on obtaining nearly maximum turbine efficiency and 
a parametric exploration of the trade-offs in turbine design, which may be 
available in the cycle-selection process. Although the analysis is not in- 
tended as a design procedure, it is general enough to demonstrate design 
trends. 

In the analysis, aerodynamic efficiency has been retained as a dependent 
variable. The correlation of reference 4 has been employed, and the effects 
of moisture on efficiency has been neglected. Efficiency reductions due to the 
presence of moisture are recognized, but adequate data are not available to 
predict their magnitude. It was assumed that the losses due to moisture will 
not be changed much by the variations investigated and that trends in "aero- 
dynamic” efficiency will reflect the trends in overall efficiency. 

An IBM 7094 computer was employed to determine the effects of geometry, 
inlet temperature, temperature ratio, stress level, power level, and the work- 
ing fluid on axial-flow turbines. 



METHOD OF ANALYSIS 


A sketch of the turbine model is 
shown in figure 1* Conceptually, it is 
a conventional multistage turbine; the 
blades of each wheel are attached to 
square rims, with the disks shaped for 
constant stress- One of the main geo- 
metrical variations was stage-to-stage 
blade mean-section diameter* Symbols are 
defined in appendix A, while details of 
Figure l. -Analytical turbine model. the analytical technique are presented in 

appendix B. In order to calculate the 

required turbine work and weight flow, the overall turbine and generator effi- 
ciencies were assumed to be 0.8 and 0.9, respectively. Efficiency reductions 
due to moisture were not evaluated, but the assumption of an overall turbine 
efficiency of 0.8 was an attempt to account for these losses. 

For the purposes of this investigation, a set of reference values has been 
chosen and is presented in table I. Each of these parameters is discussed to 
bring some of the generalizing assumptions to the reader 1 s attention. 



Working Fluid 


Arbitrarily, potassium has been employed as the basis for exploring varia- 
tions in the other turbine parameters; however, the effects of cesium, rubid- 
ium, and sodium on turbine performance characteristics have also been investi- 
gated. 

TABLE I. - REFERENCE VALUES AND RANGES OF TURBINE 


PARAMETERS INVESTIGATED 


Parameter 

Reference 

Other values 


value 

investigated 

Working fluid 

Potassium 

Cesium, rubidium, 
and sodium 

Inlet temperature, T-^, °R 

2500 

2400 and 2300 

Turbine taper ratio, Dm,n/Dm,l 

leO 

0.8 to 1.5 

Temperature ratio, Tq/T-^ 

0.772 to 0. 778 

0. 65 to 0. 85 

Allowable stress, a a , percent 

50 

40 and 60 

of rupture stress 



Inlet quality, 

1.0 


Last-stage hub-to-tip radius 



ratio, (r tL /r t ) n 

0. 7 

0.6 to 0.8 

Stator exit angle, a, deg 

70 

60 to 80 

Disk taper ratio limit. 

2.5 

2.0 to 3.0 

(^a/^d)max 



Mean blade aspect ratio 

2.0 


Blade-root stress correction 

0.7 

1.0 

factor, Tjr 



Power level, P, Mw 

1 

2, 5, and 10 

Fluid state 

"Frozen 11 

"Equilibrium” 


Inlet Temperature 

An inlet tempera- 
ture of 2500° R has 
been assumed to be a 
reasonable temperature, 
which Rankine cycle 
technology may attain. 
The effects of reducing 
the inlet temperature 
to 2400° and 2300° R 
were investigated. 


Turbine Taper Ratio 

The turbine taper 
ratio, defined as the 
ratio of the last- to 
first -stage blade mean- 
section diameter 
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®m,n/^m,l.» ^ as "been adopted as one of the primary variables. Overall turbine 
geometry was determined by satisfying the continuity equation, specifying 
(l) a fixed axial blade chord, (2) the last-stage hub-to-tip radius ratio 
(r^/ r^) n , and ( 3 ) the turbine taper ratio, and allowing the intermediate- stage 
diameters to vary linearly between the first- and last-stage diameters. Tur- 
bine performance will be compared to the reference taper ratio of 1 (constant 
mean diameter). The range of this variable included unusual turbine geom- 
etries, the last-stage mean diameter being as small as 0.8 that of the first 
stage, and somewhat more conventional variations, the last-stage mean diameter 
being as much as 1.5 that of the first stage. 


Turbine Temperature Ratio 

The reference values of turbine temperature ratio were those that yielded 
minimum specific radiator area (sq ft/kw) for the working fluids and an assumed 
condensing radiator (see appendix B and ref. 3). This parameter is a measure 
of the turbine specific work. It is important during the cycle-selection pro- 
cess because trade-offs between the size and weight of the radiator and the 
remainder of the system are heavily dependent upon this ratio. Effects of 
this ratio on the turbine configuration have been included. 


Stress Level 

The molybdenum-based alloy, containing 0.45 percent titanium and 0.07 per- 
cent zirconium (referred to as TZM), has been assumed to be the turbine mate- 
rial. It is felt that this material is typical of the refractory alloys that 
appear to be applicable. For the reference condition, 50 percent of the rup- 
ture stress over a 20,000-hour life has been used for allowable stresses. 
Effects of small variations (40 to 60 percent of rupture) were included in this 
study to indicate the sensitivity of the results to the assumed stress. 


Inlet Quality 

Saturated vapor was assumed at the turbine inlet. Reference 2 has shown 
that superheating this flow does not appear attractive because large amounts 
are required to substantially reduce moisture content while cycle efficiency 
is reduced. 


Last-Stage Hub-to-Tip Radius Ratio 

A value of 0. 7 has been assumed for the reference condition. This was 
felt to be a reasonable value on the basis of existing gas turbine technology; 
however, variations from 0. 6 to 0. 8 were also included. 
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Stator Exit Angle 

An angle of 70° has been shown to be a reasonable value for good effi- 



ciency (ref, 4) and was used for most of the calculations. This value was as- 
signed to the blade mean-section radius of the first and last stages. 
Intermediate-stage stator exit angles were determined from the axial velocity 
component, which satisfied continuity, and the tangential component, which 
satisfied a requirement of having a constant stage speed-work parameter. Vari 
ations in the assigned first- and last-stage values from 60° to 80° have been 
included in the analysis. 


Disk Taper Ratio 

A disk taper ratio of 2.5 is used here as a reference value. In general, 
this ratio should be as large as possible without causing stage spacing prob- 
lems near the axis of rotation. Limits on this ratio from 2. 0 to 3. 0 were 
investigated. Each wheel was conservatively assumed to attain the total tem- 
perature of the vapor entering that stage. Each disk was < then allowed to 
assume the value of disk taper ratio prescribed by the allowable stress and 
rotational speed. Hence, all disks were fully stressed and of minimum size. 


Mean Blade Aspect Ratio 

A single value of *2. 0 has been used. In an attempt to describe reason- 
able turbine geometries, a fixed axial blade chord was determined for each 
turbine calculation. The procedure was such that the individual stage aspect 
ratio varied below 2.0 in the front stages and above 2.0 in the last stages 
(see appendix B for more details). 


Blade-Root Stress Correction Factor 

A value of 0. 70 has been used for the reference condition. Linear vari- 
ation of the blade cross-sectional area from hub to tip has been assumed. 
Conceptually, the airfoil thickness rather than the rectangular planform was 
assumed to vary. A stress correction factor of 1 has also been investigated. 


Power Level 

An output of 1 megawatt with an assumed generator efficiency of 0. 9 has 
been used for the reference value. Overall system power levels up to 10 mega- 
watts were also considered. 


Fluid State 

For ease of analysis, the vapor fluid state has been assumed to be 
"frozen" (as used in ref. 5 with varying vapor molecular weight) for the ref- 
erence condition. Frozen properties are calculated in reference 5 as the sum 
of the contributions from the monatomic and diatomic molecules present at a 
given pressure and temperature. Not included is the energy required to alter 
the degree of equilibrium as a result of dissociation. Reference 5 also cal- 



Last-stage blade mean- Rotational speed, N, rpm Stage speed-work parameter, A Static efficiency, r\< 

section diameter, D m n , in. 




(a) Static efficiency, stage speed-work parameter, rotational speed, and 
last-stage blade mean-section diameter. 

Figure 2. -Turbine characteristics as functions of turbine taper ratio. 
(See table I for reference conditions.) 


culates "equilibrium” values 
of vapor specific heat and 
speed of sound, which include 
this energy. 

The degree of dissocia- 
tion occurring in the turbine 
expansion process is not known. 
In general, it is felt that the 
descriptive fluid state will be 
somewhere between the "frozen" 
and "equilibrium" states. 

Hence, for comparison, turbine 
characteristics using the equi- 
librium fluid state are also 
presented (see appendix B). 


RESULTS MD DISCUSSION 

The results of this para- 
metric analysis are presented 
in two sections, namely. Geo- 
metric Variations and Thermo- 
dynamic Variations. It must be 
emphasized that the magnitudes 
shown in the results are de- 
pendent upon several general- 
izing assumptions. The trends 
and relative changes will be 
stressed in the discussion. 


Geometric Variations 

Effects of these param- 
eters on turbine characteris- 
tics were investigated assum- 
ing potassium as the working 
fluid at an inlet temperature 
of 2500° R. Saturated vapor 
at the inlet and the "frozen" 
fluid state were assumed. The 
turbine temperature ratio was 
set at 0.772 (see appendix B) , 
while the generator output was 
assumed to be 1 megawatt. (See 
table I, p. 3, for a complete 
list of reference conditions. ) 


Turbine tapering. - Turbine characteristics are shown in figure 2(a) as 
a function of turbine taper ratio Ifo,nAb,l* Overall static efficiency, those 
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values of stage speed-work 
parameter that satisfy the 
turbine specific work require- 
ment, rotational speed, and 
last-stage mean diameter, are 
presented for 3- to 14-stage 
solutions. Variations of the 
last -stage rotor exit Mach 
number as calculated at the 
blade mean-section diameter, 
the maximum rotor tip speed, 
and the first-stage blade 
height are superimposed on the 
results. 

In general, for any fixed 
number of stages, efficiency, 
stage speed-work parameter, 
and rotational speed initially 
increased as the first-stage 
blade mean- section diameter 
was decreased with respect to 
the last- stage blade mean- 
section diameter. For the 
three-, four-, and five- stage 
solutions, maximum values of 
efficiency and stage speed- 
work parameter were obtained 
at values of the turbine taper 
ratio above 1. 3. Variations in last-stage mean diameter, for a fixed number 'of 
stages, were small over the range in turbine taper ratio from 0.8 to 1.5. 

Below turbine taper ratios of 1. 3, the first- stage disks of all solutions were 
stress limited. Increasing the taper ratio above 1.3, for the three-, four-, 
and five- stage solutions moved the stress limit continuously in the turbine 
model, either to the blades or to the later- stage disks. Increased tapering 
for a fixed number of stages was accompanied by increasing rotor tip speeds, 
first- stage blade heights, and rotor exit Mach numbers. The latter effect was 
small except at large values of turbine taper ratios (above 1.3). 

Figure 2(b) summarizes the effects of turbine tapering by comparing the 
peak efficiency solutions (or the maximum for n > 5, namely, those for 
A =s 1.0) to corresponding constant-mean- diameter solutions as a function of the 
number of stages. Rotational speed is divided by the n-stage constant-mean- 
diameter rotational speed Improved static efficiency or a reduction in 

the number of required stages may be obtained by permitting decreases in the 
first-stage diameter with respect to the last stage. For the present condi- 
tions, the largest improvement in efficiency was 0.06 for the five-stage solu- 
tions with the fifth- stage mean diameter 1. 4 that of the first stage. The 
corresponding rotational speed was 1. 735 that of the const ant- me an -diameter 
speed. Furthermore, It may be seen that if, say, 99 percent of the maximum 
efficiency (0. 870) were a design goal, then 11 stages would be required for a 
constant-mean-diameter solution as compared to six stages with tapering. 
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(b) Cross plot of efficiency trends. 

Figure 2. - Concluded. Turbine characteristics as functions of turbine 
taper ratio. (See table I for reference conditions. ) 


7 








Stator exit angle, a, deg 


(a) Static efficiency, stage speed-work parameter, and rotational speed. 

Figure 3. - Constant-blade-mean-section-diameter turbine characteristics as 
functions of stator exit angle. (See table I for reference conditions. ) 


The cross-hatched areas 
in figure 2(h) show a region 
of design variation in rota- 
tional speed and tapering 
available within one point of 
peak efficiency. It can be 
seen that the combined 
effects of tapering of the 
turbine flow passage and add- 
ing stages might be employed 
to afford good efficiency 
over a wide range of design 
rotational speeds. For ex- 
ample, if a static efficiency 
of 0. 80 were acceptable, 
five stages for the constant - 
mean-diameter solutions would 
be required; then if tapering 
were permissible, rotational 
speed at the design point 
could be selected anywhere 
within 1.0 and 1. 74 of the 
constant-mean-diameter speed. 
Of course, increasing or de- 
creasing the number of stages 
would permit a still wider 
range in the choice of rota- 
tional speed. 

Stator exit angle . - The 
effect of stator exit angle 
on the constant-mean-diameter 
turbine characteristics is 
shown in figure 3(a). Values 
of the stator angle that re- 
sulted in maximum static ef- 
ficiency were a function of 
the number of stages and were 
characterized by decreasing 
angle with increasing number 
of stages. For any fixed 
number of stages, increasing 
the assigned stator exit 
angles resulted in decreasing 
rotational speed, stage 
speed- work parameter, and 


last- stage rotor Mach number, while the first- stage blade height increased. 


Figure 3(b) summarizes the effects shown in figure 3(a). In general, the 
solutions for lower numbers of stages were more sensitive to stator exit angle; 
however, over the range from 3 to 12 stages, the reference value of 70° fell 
at least within 0.02 of the peak stage efficiency. Although not shown, the 
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optimum stator exit angles were also func- 
tions of the turbine taper ratio, which, 
for a fixed number of stages increased 
with increasing taper ratio. Again, the 
solutions for lower numbers of stages were 
more sensitive to this effect. 

Looking at this parameter in another 
light, it can be seen that, roughly, a 
±5° change in stator exit angle about the 
optimum value can permit a ±10 percent 
change in design rotational speed with 
little penalty in efficiency. 

Disk taper . - Effects of varying the 
disk taper ratio limit from the reference 
value are presented in figure 4, for 
constant -mean-diameter solutions. For 
any fixed number of stages, increasing 
the disk taper ratio limit resulted in 
increased efficiency and design rotational 
speed. 


(b) Cross plot of peak efficiency trends. 

Figure 3. - Concluded. Constant-blade-mean-section- 
diameter turbine characteristics as functions of stator 
exit angle. (See table I for reference conditions. ) 



Again, if 99 percent of the maximum 
static efficiency (0.865) for constant 
mean-diameter turbines were set as a goal, 
then for the reference disk taper ratio 

of 2.5, nine stages would be re- 
quired. Disk taper ratio limits 
of 3. 0 and 2. 0 then would result 
in a one- stage reduction and a 
four-stage increase in the number 
of required stages, respectively. 


Figure 4. - Effect of disk taper ratio limit on constant-biade-mean -section - 
diameter turbine characteristics. (See table I for reference conditions. ) 


Last-stage hub -to -tip radius 
ratio . - The effects of this 
parameter on constant-mean- 
diameter turbine characteristics 
are presented in figure 5. For 
any fixed number of stages, de- 
creasing the last -stage radius 
ratio resulted in increased ef- 
ficiency and rotational speed. 

If 99 percent of obtainable 
efficiency (0.867) were pre- 
scribed, then the reference last- 
stage radius ratio of 0. 7 would 
require 10 stages, while a radius 
ratio of 0. 8 would require 
14 stages and a radius ratio of 
0. 6 would require 8 stages. 
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Blade -root stress 
1 correction factor, 
* 



Figure 5. - Effect of last-stage hub-to-tip radius ratio 
on constant-blade-mean-section -diameter turbine 
characteristics. (See table I for reference condi- 
tions. ) 

Blade- root stress correction 
factor . - The effect of not varying 
the hlade cross-sectional area 
(i|r = 1. 0) is compared to the refer- 
ence condition of ^ = 0. 7 in fig- 
ure 6. As would be- expected, only 
those solutions that were near or at 
a stress limit in the blades were 
affected, namely, only the highly 
tapered, low-number-of-stage solu- 
tions (n < 7). For these solutions with = 1. 0, the stress limitation 
shifted from the disks to the blades at lower values of the turbine taper 
ratio. Hence, allowing the rotor blade cross-sectional areas to vary may per- 
mit increased design rotational speed and small increases (not more than 0.02) 
in aerodynamic efficiency for small-number-of- stage turbine designs. 
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Figure 6. - Effect of blade-root stress correction factor on turbine 
characteristics. (See table I for reference conditions.) 


Stress level . - The effect of small changes in allowable stress on turbine 
characteristics is summarized in figure 7. Variations for 40, 50 (the refer- 
ence level), and 60 percent of the rupture stress are shown. The constant- 
mean-diameter solutions for any number of stages exhibited small changes in 
design rotational speed (roughly ±5 percent) and efficiency (within ±0.02) with 
the ±10 percent change in allowable stress. 
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Turbine taper ratio, Rotational speed, N, rpro Static efficiency, 
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Figure 7. - Effect of stress level on turbine characteristics. 
(See table I for reference conditions.) 
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Figure 8. - Effect of inlet temperature on turbine char- 
acteristics. (See table I for reference conditions. ) 


The crossover among the 


The peak n-stage efficiency solutions 
were characterized by nonuniform variations 
in rotational speed with stress level for n < 7. 
three rotational speed curves was due to the combined effect of the varying 
degrees of tapering required for peak efficiency and the nonoptimum values of 
stator exit angles. Except for the apparent small improvements in obtainable 
efficiency for the three-,, four-, and five-stage solutions, the small changes 
in allowable stress did not affect the efficiency obtainable with turbine 
tapering. The stress changes did, however, alter the amount of tapering re- 
quired to obtain peak efficiency. 


Thermodynamic Variations 

Effects of these parameters on turbine characteristics were investigated 


11 









(a) Static efficiency, stage, speed-work param- 
eter, and rotational speed. 


with the reference values (see table I) of the 
geometrical restraints discussed previously. 

Inlet temperature . - Characteristics of 
the turbine solutions over the range in inlet 
temperature from 2500° to 2300° R are summa- 
rized in figure 8. Comparison of the effi- 
ciency variations shows that, for a fixed num- 
ber of stages, the improvement from constant- 
mean-diameter solutions to tapered turbine de- 
signs decreased with decreasing temperature. 

The degree of tapering (except for n < 5) and 
relative change in rotational speed required 
to achieve peak efficiency, however, decreased 
with decreasing inlet temperature. Again, the 
low-number-of-stage solutions exhibit less 
uniform variations because of their greater 
sensitivity to the assigned stator exit angle. 

Variations in maximum rotor tip speed are 
included in this figure to give an indication 
of the relative increases due to tapering. 

Over this range in inlet temperature, there 
appears to be a narrow band of tip speeds that 
is a function of the number of turbine stages, 
which, if attained by a design, will poten- 
tially yield nearly peak efficiency. 

A 0. 01 decrease occurred in the maximum 
efficiency obtained between inlet temperatures 
of 2500° and 2300° R. This was largely due to 
a corresponding decrease in Reynolds number; 
however, the relative effect of inlet tempera- 
ture may be compared by again setting an arbi- 
trary goal of 99 percent of the maximum effi- 
ciency at each temperature. Then, at 2500°, 
2400°, and 2300° R, the reduction in the number 
of required stages from the constant-mean- 
diameter to tapered solutions was four, three, 
and two, respectively. Hence, turbine tapering 
was more effective in reducing the number of 
stages required for nearly peak efficiency as 
inlet temperature increased. 


Figure 9. - Constant-blade-mean -section- 

diameter turbine characteristics as functions Turb ine temperature ratio . - Effects of 

of turbine temperature ratio. (See table I this variable on turbine characteristics are 

for reference conditions.) . * \ 

shown in figure 9(a) for constant-mean-diameter 
solutions. Increasing turbine temperature 
ratio from 0. 65 to 0. 85 corresponds to decreasing the specific work by 60 per- 
cent and increasing the exit quality by 14. 6 percent and required weight flow 
by 150 percent (see table II, p. 22 ). For a fixed number of stages, the aero- 
dynamic efficiency (i.e. , static) increased with increasing temperature ratio. 
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Rotational speed, N, rpm Static efficiency, 7 ? 
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Since the efficiency reduc- 
tion due to moisture could be 
expected to decrease as tur- 
bine temperature ratio is 
increased, overall turbine 
efficiency would increase 
faster than indicated in 
figure 9(a). 

For a fixed number of 
stages, increasing the tur- 
bine temperature ratio also 
increased the stage speed- 
work parameter and the first- 
stage blade height, while the 
rotational speed increased, 
reached a maximum, and then 
decreased. (Although not 
shown, the mean diameter of 
the solutions minimized where 
n was maximized. ) The tem- 
perature ratio at which the 
rotational speed maximized 
was a function of the number 
of stages and increased as 
the number of stages de- 
creased. 

Changes in efficiency and rotational speed are compared in the cross plot 
figure 9(b). The obtainable maximum static efficiency (neglecting moisture 
effects) decreased with decreasing turbine temperature ratio. Also, consider- 
ably more stages were required to obtain a given level of efficiency as turbine 
temperature ratio was decreased; however, turbine tapering could be employed 
to reduce the number of required stages. 

Selection of low turbine temperature ratios may penalize the system with 
reduced turbine performance and the complexity of many turbine stages. 

Power level. - Since a change in generator power level is reflected in the 
turbine by a linear change in weight flow, it is reasonable to expect that tur- 
bine diameter and rotational speed should scale with a change in power level. 
That is, for any given turbine taper ratio and number of stages, the diameters 
should vary directly as the square root of the power level, while the rota- 
tional speed should vary inversely as the square root of the power level. 

Turbine solutions were obtained over the range in generator power level 
from 1 to 10 megawatts at the reference conditions of table I. Within a 
±2 percent variation and independent of number of stages and turbine taper 
ratio, the results correlated with the scaling procedure. Thus, for example, 
if a turbine design for a 1 megawatt generator output prescribed a last-stage 
mean diameter of 1 foot and a rotational speed of 20,000 rpm, then the similar 
turbine configuration at 10 megawatts of generator output would have a diameter 




(b) Cross plot of efficiency and rotational speed trends. 

Figure 9. - Concluded. Constant-blade-mean-section-diameter turbine charac- 
teristics as functions of turbine temperature ratio. (See table I for reference 
conditions.) 
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Figure 10. - Characteristics of cesium tur- 
bines. (See table I for reference conditions. ) 

of 3.16 feet (+0.06), while the ro- 
tational speed would he 6320 rpm 
(±126). 

Working fluid . - Effects of 
taper ratio on turbine characteris- 
tics when cesium, rubidium, and so- 
dium are used as working fluids are 
shown in figures 10 to 12, respec- 
tively. 
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Figure 11. - Characteristics of rubidium turbines. (See table I 
for reference conditions.) 
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Figure 12. - Characteristics of sodium turbines. (See table I for refer- 
ence conditions.) 


The cesium solutions were 
characterized by the fewest 
number of stage solutions and 
very limited flexibility with 
respect to design variations. 
Three stages of nearly constant 
tip diameter (bm^n/^m^l = 0.95) 
represented the highest obtain- 
able efficiency with a reason- 
able turbine shape. This maxi- 
mum efficiency solution was 
also characterized by a fairly 
reasonable first- stage blade 
height (about 0.6 in.) and the 
lowest level of rotor tip speed 
(about 640 ft/sec maximum) 
among the peak efficiency solu- 
tions of the four working 
fluids. 

The rubidium solutions 
(fig. 11) were characterized by 
a few more number-of-stage 
solutions and a little more de- 
sign flexibility than those for 
cesium. The constant-tip- and 
constant-mean-diameter solu- 
tions for any number of stages 
were near peak efficiency. 
Either six constant-tip- 
diameter stages or five nearly 
constant-mean-diameter stages 
(Djn^n/^m^l = 0.985) afforded 
maximum efficiency. 

The sodium solutions 
(fig. 12) were characterized 
by four to twenty or more 
stages with large regions of 
design control with taper 
ratio. Up to the maximum taper 
ratio investigated (1.5), only 
the four- and five- stage solu- 
tions were approaching the 
shift in stress limitation from 
the first-stage disk. It may 
also be noted that the maximum 
efficiency level is 0. 85 in 
figure 12. In general, it is 
felt that, if more stages 
and/or larger turbine taper 
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-'ratios were investigated, the obtainable 
static efficiency level would approach 
those obtained with the other working 
fluids (approximately 0.88). 

In general, the number of stages re- 
quired for peak efficiency decreased with 
increasing molecular weight of the working 
fluid (i. e. , from sodium to potassium to 
rubidium to cesium), while the effective- 
ness of the turbine flow passage tapering 
in improving aerodynamic efficiency de- 
creased. At a turbine inlet temperature 
of 2500° R, with sodium as the working 
fluid, 12 stages tapered by at least 
50 percent would be required to attain 
maximum efficiency; while for potassium, 

7 stages tapered by roughly 30 percent 
would be required, and 5 and 3, untapered 
stages would be required for rubidium and 
cesium, respectively. 

On the other hand, it may also be ad- 
vantageous to have the larger degree of 
freedom in choice of design rotational 
speed that would be available by the com- 
bined effect of tapering and stator exit 
angle previously exhibited for potassium 
turbine solutions. 

Fluid state. - The effect of the 
"frozen" fluid ^state assumption of this 
analysis on Mach number level and rotational speed is compared with an "equi- 
librium" state in figure 13. Both constant-mean-diameter and tapered peak- 
efficiency configurations are compared. It can be seen that the frozen assump- 
tion resulted in lower Mach number levels (by about 0.1) and higher rotational 
speeds (by not more than 18 percent) as compared to the equilibrium results. 



Figure 13. - Comparison of fluid state assumptions on 
turbine characteristics. (See table I for reference 
conditions.) 


CONCLUSIONS 

Axial-flow, alkali-metal vapor turbines have been investigated as one of 
the components in megawatt space electric powerplants. A parametric analysis 
has been conducted to determine the interactions of several design variables 
in obtaining maximum turbine efficiency. Efficiency reductions due to the 
presence of moisture were not included* The following conclusions were 
reached: 

1. Turbine tapering (i.e. , changing the ratio of last- to first-stage 
blade mean-section diameter) was an effective means for attaining maximum effi- 
ciency with a minimum number of stages. Turbine taper ratios for maximum effi- 
ciency decreased with increasing working fluid molecular weight. At an inlet 
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temperature of 2500° R, these solutions were characterized hy the following: 

(a) 12 stages for sodium,, with a taper ratio of at least 1.5 

(b) 7 stages for potassium, with a taper ratio of roughly 1.3 

(c) 5 stages for rubidium, with a taper ratio of 0.985 (nearly constant 

mean diameter) 

(d) 3 stages for cesium, with a taper ratio of 0.95 (nearly constant tip 

diameter) 

2. Turbine tapering and/or small variations in stator exit angle permit a 
wider selection of design rotational speed, while still maintaining a desired 
level of efficiency. 

3. Advantages of turbine tapering in reducing the number of stages re- 
quired for nearly maximum efficiency decreased with decreased inlet tempera- 
ture. 


4. Selection of a low ratio of outlet to inlet temperature may penalize 
turbine performance. 


Lewis Research Center 

Rational Aeronautics and Space Administration 
Cleveland, Ohio, August 20, 1964 
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APPENDIX A 


SYMBOLS 

A annular flow area, sq in. 

a speed of sound, ft/sec 

b blade height, in. 

Cp vapor specific heat at constant pressure, Btu/(lb)(°R) 

c axial blade chord, in. 

D diameter, in. 

D m , n /D m ,i turbine taper ratio 

g acceleration due to gravity, 52.174 ft/sec 2 

h enthalpy, Btu/lb 

J mechanical equivalent of heat, 778.26 (ft)(lb)/Btu 

M Mach number 

m vapor molecular weight 

N rotational speed, rpm 

N/Mjj ratio of rotational speed to n-stage constant-mean-diameter rota- 

tional speed 

n number of stages 

P power level, Mw 

He Reynolds number 

r radius, in. 

r^/r-t hub-to-tip radius ratio 

S entropy, Btu/(lb)(°R) 

T temperature, °R 

t^ thickness of disk at rim, in. 

U wheel speed, ft/sec 
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V 

absolute velocity, ft/sec 


w 

velocity relative to rotor blade, ft/sec 


w 

weight flow, Ib/sec 


X 

quality of vapor 


y a 

thickness of disk at axis of rotation, in. 


y a Ad 

disk taper ratio 


a 

stator exit angle at blade mean- section radius, measured from 
direction, deg 

axial 

3 

rotor angle at blade mean-section radius, measured from axial 
tion, deg 

direc- 

T 

ratio of specific heats 


A 

difference 


n 

efficiency 


X 

stage speed-work parameter, Um 2 /gJ(Ah) 



vapor viscosity into turbine, lb/(ft)(sec) 


p 

vapor density, lb/cu ft 


Ptzm 

density of turbine material, 639 lb/cu ft 


a 

stress, psi 



stress correction factor for tapered blades 


Subscripts : 


a 

allowable 


b 

blade 


d 

disk 


e 

equilibrium 


g 

generator 


h 

hub 

« 

i 

inlet 

♦ 
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id ideal 

j dummy variable referring to stage number, 1 to n 

k constant-blade mean-section diameter 

l liquid 

m blade mean section 

max maximum 

n last stage 

o outlet 

p peak 

s static 

T turbine 

t tip 

v vapor 

x axial component 

9 tangential component 

Superscript: 


overall 



APPENDIX B 


ANALYTICAL PROCEDURE 

The investigation was programed for computation on an IBM 7094 digital 
machine. The working equations and the general approach are presented in this 
appendix. 


Input Conditions 

Input to the program included working fluid properties, assigned turbine 
requirements, and a series of constants. All working fluid properties with 
the exception of vapor viscosity were obtained from reference 5. The proper- 
ties required for the turbine computation procedure were vapor density p, 
speed of sound a, and the ratio of specific heats y as .functions of tempera- 
ture and single values of inlet vapor viscosity |j. (obtained from ref. 6). 

Both "frozen” and "equilibrium” values (as defined in ref. 5) of the speed of 
sound and the ratio of specific heats were investigated. The equilibrium 
ratios of specific heats were calculated from the properties listed in refer- 
ence 5 as 


TA - 


-'p A e 


tp,e 


1. 9865 
nu 


(Bl) 


Interpolation of a, y, and p in the machine program was done by the Gregory- 
Newton method (ref. 7, p. 264) by using three forward differences. Because of 
the large rate of change of density with temperature, this property was Inter- 
polated logarithmically by using the Gregory-Newton method. 

The assigned turbine requirements were inlet temperature T-^, temperature 
ratio T 0 /Tj_, overall specific work Ah, weight flow w, and exit quality Xq. 
For each fluid and inlet temperature, those turbine temperature ratios that 
yield minimum specific radiator area (sq ft/kw) were determined from the method 
of reference 3 (modified Rankine cycle), but using an assumed generator effi- 
ciency of 0. 9. In order to calculate the required turbine work and weight 
flow, the overall turbine efficiency was assumed to be 0.8. Then, the overall 
turbine work was evaluated from 

Ahrp = ^(hj, - ho^idh 

= 0.8(h 1 - h Q ^ ld ) (B2) 

where 

ho, id ~ il o,v ” ^(So^ - Sj_^ v ) (B3) 

the turbine weight flow from 
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w = 949 — — — = 1053 =J- 

rigAhij Ahiji 

and the turbine exit quality from ■ 

h o - h o,Z 

k<D,V " ^ 0,7 

For reference, these turbine working conditions are given in table II. 


TABLE II. - TDKBME WORKING COEDITIONS 


Inlet 

Work- 

Temper- 

Overall 

Outlet 

Specific 

Inlet 

temper- 

ing 

ature 

turbine 

quality. 

weight 

viscosity. 

ature. 

fluid 

ratio. 

work. 

x 0 

flow. 




To/Ti 

Ahjji, 


lb/(sec) (Mw) 

lb/ (ft) (sec) 

°R 


Btu/lb 




2500 

Sodium 

0. 772 

302 

0. 854 

3.49 

1. 53X10“® 


Potassium 

.776 

150 

.877 

7.02 

1. 70 


Rubidium 

.778 

60.0 

.889 

17.5 

2. 03 


Cesium 

.778 

36.0 

.896 

29.3 

2.16 

2500 

Potassium 

0.65 

245.5 

0.801 

4.29 

1. 70X10" 5 



. 70 

207.0 

.831 

5.08 




! .75 

169.4 

! .861 

6.21 




.80 

132.9 

.891 

7. 92 




.85 

98.2 

.919 

10. 71 


2400 

Potassium 

0. 775 

152 

0.872 

6.92 

1.50X10" 5 

2300 

Potassium 

.775 

154 

.869 

6.83 

1. 47X10-5 


Other pertinent input quantities included the assumed turbine material 
density PipzM* allowable stresses u a , blade mean-section stator exit angle a, 
last-stage hub-to-tip radius ratio (rj 1 /r^) ri , turbine taper ratio E^ h/I^ q, 
and disk taper ratio limit (y a Aq) max . Arbitrarily, a molybdenum-based alloy 
containing 0.45 percent titanium and 0.07 percent zirconium (TZM) has been 
assumed as the disk and blade material. Allowable stresses were determined 
from the data of reference 8 by extrapolating rupture stress on a Larson-Miller 
plot to 20,000 hours. In general, 50 percent of these extrapolated rupture 
stresses were used as the allowable stress. These were 0 a = 9. 8X10^ pounds 
per square inch at 1600° R and a a = 9. 9X10^ pounds per square inch at 2700° R 
with an assumed linear variation of the logarithm to the base 10 of the 
stresses between these temperatures. 


(B4) 


(B5) 


Work- Speed Balance 

The first portion of the calculation is to determine all combinations of 
n. A, and N in the range 0. 2 < A < 1. 0 and 2 < n < 20 that satisfy the 
overall turbine work requirement within the assumed allowable stresses. A 
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solution to a problem is initiated by setting A = 0. 2 and n = 2. An itera- 
tive procedure is then employed to find the compatible values of the axial 
components of velocity into the first- and last- stage rotors (V x i 1 and 
Vx,i,n* see fig. 14 for velocity diagram nomenclature) and rotational speed 
for the minimum value of the stage speed-work parameter. 



Figure 14. - Velocity diagram nomen- 
clature. 


Initially, arbitrary values are assigned to 
V x i i and V x ^q^ n , and the temperature and quality 
of the working fluid into the last stage are estimated. 
Then, from the geometry, the absolute velocity into 
the first and last stage at the mean radius is given by 

v i = v x,i V 1 + tan 2 a- ( B6 ) 

Then, in order to correct to static flow conditions, 
the Mach numbers and the density Ratios (assuming 
perfect gas relations) are evaluated from 



where a is the speed of sound corresponding to the total temperature Tj_ 


entering the stage in question, and 

1 


(?) - ( 


(B8) 

\ w s/ ± 




The annular flow area required to pass only the vapor is calculated from 


Since (r h /r t ) n 
follow from 


A = 


144V. /p\ 

i V x,i V>J ± 


(B9) 


is specified as input to the problem, the last- stage radii 



(BIO) 


r m,n 


r 


t ,n 


J 
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The first-stage radii follow from the prescribed taper ratio Dm nA(m,l and 


e , _ 

Pm,n 

%,1 


~\ 


r tA 

r h,l 


r m,l 

%,1 


+ 4ltr m,l 
Al 

" 4irr m,l 


> 


J 


(Bll) 


Intermediate-stage radii are calculated proportionately between the first- and 
last-stage radii. In order to determine the radii of the disks at the rims, a 
mean axial chord is determined from 


where 


c 


c, + c 
1 n 


2 


(B12) 



_ r t,m ~ r h,n _ ^n 
2 2 


with the arbitrary restriction that c^ and c n be equal to or greater than 

l/2 inch. Then, since square rims were assumed (of width equal to c), the 
radii of the disks at the rims are calculated from 


r 




(B13) 


where j = 1, 2, . . . , n 


Angular speed is now determined as the minimum value from among those 
calculated by assuming the blades and then the disks of each stage to be stress 
limited. The radial stress at the hub of tapered blades is given by (see 
ref. 9) 


2(l2) 4 (30) 2 g 


(B14) 


Therefore, if the blades are stress limited, 

N = 4320 


4 


2go, 


a 


^TZM 71 


(B15) 
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Conservatively, it was assumed that the steady-state temperature of >the wheels 
was uniform and equal to the temperature of the vapor entering the stage. 
Hence, a a is evaluated for each stage at T-j^j. The stress level for a 
constant-stress disk with no temperature gradient is evaluated from (see 
ref. 9) 


°d 


jt 2 p TZM N' 2 r| 


2(l2) 4 (30) 2 g In 



(B16) 


Therefore 



, if a disk is stress limited, it attains the disk taper ratio limit 
assigned to the problem, and the rotational speed is 


N = 



(B17) 


The next step in the procedure is to check the resulting first- and last- stage 
stator exit angles at the mean radius (cc^ and a^) against those assigned by 
the problem. The stage speed- work parameter is defined (ref. 4) as 


and may also be written as 

A = 


gJ Ahj 


u m . rtNr m 


AV 0 , m “ SOCVe^i - V 0 , m , o ) 


(B18a) 


(B18b) 


Since A is specified and N and r m have been calculated, the change in 
the tangential velocity component for the first and last stage are computed 
from equation (B18b). Also, depending on the current value of A, Vq ? ± is 
found from 


V, 


M = (* + l) AV « 


for 0 < A < 0. 5 -1 


Vq } ± = AVg for 0.5 < A < 1.0 


(B19) 


Then, 


tan a = 


^8, m, i 


V- 




(B20) 


If either or both the first- and last- stage stator exit angles do not agree 
with the assigned values, then new axial components of the absolute velocity 
are computed as 


I 
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(B21) 


V. 


x,i 


v e,m,i 
tan a 


Before Being iterated., the stage and total turbine work are calculated 

from 


U . 

ZNh, = JSjA 
J gJA 


3 = 1, 2, 


n 

Ah^ = ^ Ahj 
0=1 





(B22) 


and new values are computed for and Xq^ n . The change in temperature 

and quality across a turbine stage was assumed to be proportional to the stage 
work. Therefore, 


T- 


i,n 


= T, 


o.n 


X i, 


= X, 


■o.n 



Ah n 



n 

ATip 


Ahiji 


I 

> 


Ab n 


+ 

AXjrp 


Ahf 



(B24) 


The procedure now goes hack to equation (336) and iterates until the cor- 
rect values of V x anc ^- n are obtained for the assumed A = 0.2. 

Now if the total turbine work is greater than the required turbine work, a new 
value of A is calculated as the current work divided by the required work 
times the previous value of A. The procedure again reverts to equation (B6) 
and iterates until the turbine work equals the requirement. On the other hand, 
if the current turbine work for A = 0. 2 is less than the required amount, 
then no solution is possible and the number of stages is increased by one. The 
calculation then repeats by assuming arbitrary values of Vx,i,l and v x,i,n 
and with the new value for n goes back to equation (B6). 

Once the work and speed requirements are balanced as described in this 
section, the calculation continues as in the following section. After the 
entire calculation for a particular n is completed, a new calculation is 
made for n + 1 . A problem is terminated either when A > 1 . 0 or n > 20. 


Efficiency Estimation 

The approach of Stewart (ref. 4) was used but modified to include the 
relatively small effect of reheat (enthalpy level). The entire procedure is 
included for the reader T s convenience. Stage total and static efficiencies 
are given, respectively, in reference 4 as 
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(B25) 


Ah _ A 
^id A + i A 


and 

T, - ^h - h 

S ^id.s A + |(A + B) 

where 

A « 0,4Re ill (C + 2D) 
cot a 

with the Reynolds number given by 


Re = 


12w 


and the expressions B, C, D evaluated in table III. 


(B26) 


(B27) 


(B28) 


TABLE III. - VARIATION OF EFFICIENCY TERMS 


WITH STAGE SPEED- WORK PARAMETER 



Overall static efficiency may be expressed as 


_ £ Ahj 

tJ = _ Ah — (B29) 

S ^kldjs ^id^s 

Also, 
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n-1 

^id^s = + ^id^s^n + “ ®i,l) (B30) 

J-l 


now. 


Ah 


id,s,n 


1 

%,n 


Ah n 


(B31) 


and 


with 


n=l 

Sl,n - S 1A = £ (S DjJ - S 1>} ) (B32) 

j=l 


(Ah± d - Ahjj - T 0 ^j(S 0 ^ t 


“ B i,J> 


(B33) 


or 


n-1 



(B34) 


Combining equations (B29) to (B3l) and (B34) yields the overall turbine static 
efficiency as 


n s - 


E 

J-l 


Ah, 


s - n-1 n-1 

E + ^ Ah ° + T °’ s > n E — (i- nAh, 


(B35) 


J-l 


33. T °'J V'J 


The value of the static exit temperature is calculated from the ideal gas 
relation: 


T 0 ,s,n - To,n(l + ^ <n) 


-1 


(B36) 
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